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Abstract 

We calculated ah initio the three-fold differential cross section of a double single-photon Helium 
photoionization at the equal energy sharing, and obtained from one the Gaussian width parameter 
7, describing the angular interelectron correlations, for the total electrons energies range E from 
0.1 eV to 100 eV. The results are in the excellent agreement with experimental data but indicate 
that the Wannier threshold law for the angular distribution 7 oc 

£1/4 

is not correct at energies 

attainable in modern experiments. It is shown that the 7 dependence on the energy is much better 
described by the modified threshold law, obtained by Kazansky and Ostrovsky [J. Phys. B 26, 
2231 (1993)]. Also, we explored the Gaussian width parameter for a double photoionization of the 
targets with the strongly asymmetrical initial state configuration: the atomic Hydrogen negative 
ion H~ and the Helium in the 2s 1 S and 3s 1 S excited states. We found that the Gaussian width 
dependence on the total ejected electrons energy for these targets has a maximum at low energies. 
We show also that the correlation parameter dependence on the interelectron angle for these targets 
is essentially non-Gaussian and has a number of peaks equal to a number of initial state radial 
nodes, that reveals the new abilities for the qualitative analysis of the electron structure. 

PACS numbers: 32.80.Fb,32.80.Gc 
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I. INTRODUCTION 



It is known, that if the incident radiation is linearly polarized in the Oz direction, the 
three-fold differential cross section (3DCS) of a double photoionization of an atom by a 
single photon can be represented via gerade and ungerade amplitudes []J 
d 3 a 

, F , , = £ 2 ,0i 2 )(cos0i + cos# 2 ) + a u (E 1 ,E 2 ,9 12 )(cos9 1 -cosfl 2 )| 2 , (1) 

dEid\lid\l 2 

where the ungerade amplitude a u = for electrons energies E\ = E 2 , and the gerade ampli- 
tude usually called a correlation parameter can be approximated at low electrons energies 
by a Gaussian curve 



a g (E 1 ,E 2 ,6 12 ) ~ Aexp 



-21n2- 



ry2 



(2) 



as it was shown by Rau [2] following Wanniers theory |3J. The Gaussian width parameter 7 
is a single angle parameter describing the angular distribution, therefore it is often used for 
the analysis of strength of the angular interelectron correlation. When 7 is large then the 
interelectron correlation is weak and vice versa. According to the Wannier threshold law, 
near the double photoionization threshold it should be 

7 = 7o£ 1/4 (3) 

where E = E\ + E 2 is a total ejected electrons energy. In spite of the fact that the energy 
range of the correctness of this law is not established yet, experimentators and theorists still 
use the expression ([3]) for the data interpretation, trying to find the scaled width parameter 
7o- A lot of formulas for 70 have been proposed by various authors 4]. However, Kazansky 
and Ostrovsky have shown [5], that the Wannier threshold law (J3j) is correct only when the 
electrons deceleration by the nucleus field is neglected or for extremely law energies of the 



order of eV. From the other hand, we has found in [6] that 7 for the Hydrogen negative 
ion H - starts to grow with energy decreasing at low energies in the obvious contradiction 
with the Wannier threshold law. The aim of the present work is a calculation of the Gaussian 
width 7 at the smaller energies. 

II. DETAILS OF THE CALCULATION PROCEDURE 

In our calculations we used the time-dependent scaling (TDS) method . The main 
advantage of this method is the 3DCS obtaining for all ingoing photon energy values by 



the one act of computation. The problem is the presence of not only the wave function 
component corresponding to the double ionization and being described correctly at the 
extending coordinate system, but also components describing the bound states and the 
single ionization states which are described incorrectly for the large values of the expansion 
coefficient a(t) (see [7|). Therefore at evolution times t ^> 1000 these states generate the 
noise looking like short-period oscillations with the wavelength of the order of the radial 
grid step h. This problem requires the circumspective approach because some oscillations 
in 'y(E) appear also in calculations jsj]. But all indicates that oscillations in our results are 
the numerical artifact and come from the bound states destruction: its wavelenght is always 
of the order of the grid step for an any grid step value choice, ones appear when the a(t)h 
becomes of the order of the typical bound state radius and spread toward to the increasing 
radius. For this reason we filtered the wave function after the evolution by an eliminating 
all components with wavelenghts less than Ah. We should note that there are no oscillations 
in our calculations even without the filtration at the energy range where oscillations in C 
calculations [8] appear. 

In the calculations of the Helium photoionization given below, we used the numerical 
scheme parameters (see [6| for details) as follows: an angular basis parameter /2max=13, a 
uniform radial grid having a number of the radial nodes iV r =500 and a size £ max =25., a 
complex scaling radius £ sc =22.5, a complex scaling angle # sc =30°, a grid expanding rate 
0oo=O.l, the evolution was simulated up to the time t max = 12800. For other targets other 
radial grid parameters are used: N r =500, £ m ax=50., £ sc =45., doo=0.05 for the H~; iV r =1000, 
£max=50., £ sc =40., (200=0.05 for the He in the excited ls2s X S state; and N r =14Q0, £ m ax=70., 
£ sc =60., doo=l/30 for the He in the ls3s X S state. 

After the 3DCS has been calculated, it is necessary to obtain the Gaussian width param- 
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eter 7 from the 3DCS. The squared module of the correlation parameter \a g (Ei, Ex, 9 12 ) 
may be expressed from the 3DCS through the Eq.flT]) and, after that, fitted by (j2J) using the 
least-squares method (LS) as in the theoretical work [9|]. An alternative approach is based 
on the fitting of the twofold differential cross section (2DCS) a^ 2 \Ex, E 2 , 6x2) = dE~le~ 
the formula 



aW(Ex,Ex,e 12 ) ~ ^|A| 2 exp 



41n2(7r - 9 12 ) 2 
T 5 



2 ( A\ 

cos — (4) 



deduced from the Eq.([T]) by integrating over all angles except the #12. The analogous ap- 
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proach was used in the experimental work (lOj. Since the correlation parameter a g {9\2) can 
noticeably deviate from the Gaussian shape, the 7s calculated by two methods mentioned 
above are different. We will denote the Gaussian width parameter obtained by the fitting 
of \a g (9i 2 )\ 2 as 7(|a 9 | 2 ), and 7 obtained by the fitting of <7^(9i 2 ) as ^(a^). 

III. THE PHOTOINIZATION OF THE HELIUM IN THE GROUND STATE 
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FIG. 1: The Gaussian width parameter 7 as a function of the full energy of ejected electrons E 
for the photoionization of He: the TDS with the Gaussian fitting of |a 9 | 2 (thick solid line) and 
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In the Fig{T] we show the Gaussian width 7 as a function of the full energy of ejected 
electrons E for the photoionization of Helium in the ground state. It is obvious that our data 
coincides very well with the experiment in the whole range from 0.1 eV to 100 eV, except the 
point at 4 eV from [15(. The exact coincidence of our curve 7(0"^) with the experimental 
points at 0.116 and 0.209 eV from [10] is the most remarkable. The curve 7(|a 9 | 2 ) coincides 
with these points in rather less degree because 7(cr( 2 )) has been got in the experiment 10]. 
Generally, the differences between 7(|a 9 | 2 ) and 'j(a^) may be the feature of the degree of 
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the a g {6\2) deviation from the Gaussian shape, although the 7(|a 9 | 2 ) coincidence with the 
7(cr( 2 )) does not mean that a g (8i2) is exactly Gaussian. The Fig{T]is given in the logarythmic 
scales at both axises, in which the power dependences like the Wannier law ([3]) should look 
like sloped straight lines. Indeed, we see that our plots are close to lines when E is less than 
few eV. But the exponent is not equal to 1/4 at all! The approximation of 7(0"^) curve at 
the E range from 0.1 to 2 eV using the power law of the general form 



7 = loE s 



(5) 



through the least-squares approach yields the exponent s=0.097 and the proportionality 
constant 7o=70.2° eV _s . Such a significant deviance from the Wannier threshold law often 
used for the interpretation of experimental and theoretical data seems to be discouraging. 
However, Kazansky and Ostrovsky [5] show that the Wannier threshold law strongly modifies 
when the electrons deceleration by the nucleus field is taken into consideration because the 
electrons velocity is non-equal to zero even at zero total energy. In the Figj2]we show the 
scaled width parameter jo(E) = , y(E)/E 1 ^ in compare with the approximate curve obtained 
in j^j, which we called the Kasansky-Ostrovsky threshold law. 
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FIG. 2: The scaled width parameter 70 as a function of E for He: the TDS with the Gaussian 
fitting of I a g 1 2 (thick solid line) and (thin solid line), and the Kasansky-Ostrovsky threshold 
law 5] (dashed line). 
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FIG. 3: The Gaussian width 7 as a function of E for the photoionization of H~: the TDS with the 
Gaussian fitting of |a g | 2 (thick solid line) and (thin solid line), and the Kasansky-Ostrovsky 
threshold law [5| (dashed line). 



IV. THE PHOTOIONIZATION OF THE TARGETS WITH THE STRONGLY 
ASYMMETRICAL INITIAL STATE CONFIGURATION 

In our previous work j^l we did a comparison of ^(E) for various Helium-like ions and 
found that for the negative Hydrogen ion H~ it starts to increase at energies below 2.5 eV 
unlike other considered targets. Here we obtained the results at the energies down to 0.06 
eV. It is clear from the Figj3]that 7(cr < - 2 - ) ) increases at the energy decreasing in the range from 
2.6 eV to 0.23 eV. At energies below 0.09 eV ^(a^) is a power function of the energy with 
the exponent s=0.083 and 7o=74° eV~ s , though this range is too small for such rigorous 
conclusions. It is obvious from the Figj3] that there is a strong distinction of our results 
from the Kasansky-Ostrovsky threshold law for the nuclear charge Z = 1 p], unlike the 
Helium case (Figj2]). We should note that the "f(E) dependence obtained by Kasansky and 
Ostrovsky is monotonous (Figj3]) despite the 70 (E) dependence oscillating. Our hypothesis 
is that such distinction of the results for H~ from the results for the Helium and other 
Helium-like ions p] comes from the fact that the H~ bound state configuration is strongly 
different from the one of the Helium ground state. Indeed, when n ( 2 — > 00 the H~ bound 
state wave function has an asymptotic 
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FIG. 4: The scaled width parameter 70 as a function of E for H~: the TDS with the Gaussian 
fitting of \a g \ 2 (thick solid line) and (thin solid line), and the Kasansky-Ostrovsky threshold 
law 5] (dashed line). 

The H~ is a deuteron-like weakly bound system consisting from the Hydrogen atom and 
the electron located the most time outside the region where the attracting potential acts. 
We performed here the calculations for the another targets with the strongly asymmetrical 
initial state configuration: Helium atoms in the excited states 2s X S and 3s X S. 
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FIG. 5: The Gaussian width 7 as a function of E for the photoionization of He(2s 1 S): the TDS 
with the Gaussian fitting of |a 9 | 2 (thick solid line) and (thin solid line), and CCC results 



In the Figj5] we show 7 as a function of E for the photoionization of Helium in 2s X S 
metastable state. Our results do not deviate much from the CCC results {s| by the magnitude 
but deviate strongly by the behavior because there are no oscillations in our results. When 
the energy decreases 7 decreases at first and then begin to increase (for 'j(a^) it happens 
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at E=2.5 eV) as well as for H . Surprisingly, at E^5 eV the difference between 7(|a 9 | 2 ) 
and 7(cr( 2 - ) ) became enormous. The origin for this is clear from the Figj6j The correlation 
coefficient has a strongly non-Gaussian shape with the two peaks even at low energy E=l 
eV. When the energy decreases the secondary peak declines and therefore the distribution 
turns to Gaussian as follows from the Wannier's theory. When the energy increases the 
secondary peak grows, and becomes more than the primary one at £7=11.3 eV as it is seen 
from the Figj6b,c. Then the Gaussian fitting becomes non-applicable of course. 

In the FigJ7]we show 7 as a function of E for the photoionization of the Helium in the 3s X S 
state. Our results are strongly different from the CCC results [8] both by the magnitude and 
by the behavior. The general curve shape is similar to the one for 2s 1 S, but 7(cj( 2 )) reaches 
the local maximum at E—1.5 eV. It is seen from the FigJS]that the correlation parameter is 
strongly non-Gaussian as well as for 2s 1 S, but has three peaks at low energies. When the 
energy increases two lesser peaks grow and merge, and the cr( 2 )(# 12 ) view becomes rather 



complicated at large energies. One can see that the number of considered targets a 



(2), 



'12 J 



peaks at low energies is equal to the number of the peaks of the target "outer" electron 
density dependence on the radius in the initial state. Particularly, it is equal to one for 
H~, two for He(2s 1 S) and three for He(3s 1 S). It seems not to be an accidental coincidence 
though we still can not propose an exact explanation of this effect. If this dependence will 
be confirmed it will reveal the new abilities for the qualitative experimenatal analysis of the 
target electron structure using the over-threshold double photoionization by analogy with 
the (e, 2e) spectroscopy [18 |. 

When we have already established that the local minimum appearance in the ^{E) is 
common for all considered targets with the asymmetrical initial state configuration we can 
explore the origin of this effect. For all targets the local minima of the 7 are observed at 
the energies not depending direct proportionally on the target first ionization potential Ii, 
but of the order of 1%. We should note that despite of the absence of the local minimum 
on the plot for the Helium in the ground state (FigJTJ) there is a curve bend at E=20 eV 
approximately that indicates to some change taking place at the energy of the order of 1%. 
In the non-sequential double ionization considered here one of the electrons is ionized at 
first by the photon impact, and the second electron might be ejected through the a sudden 
change of the atomic potential (so called the shake-off mechanism), or through the first 
electron impact (so called the final state scattering). If the initial state configuration is 
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FIG. 6: The 2DCS as a function of the interelectron angle 6\2 for the photoionization of He(2s 1 S) 
for a) E=l eV; b) E=ll eV; and c) E=20 eV: "exact" TDS results (solid line), Gaussian fitting 
of cr( 2 ) (dashed line) and the Gaussian fitting of |a 9 | 2 (dotted line). 
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FIG. 7: The Gaussian width 7 as a function of E for the photoionization of He(3s 1 S): the TDS 



strongly asymmetrical then the "inner" electron momentum density is much broader than 
the one for the "outer" electron, and the single ionization cross section dependence on the 
ejection energy for the "inner" electron decreases much slower than for the "outer" one. 
It means that the process when the "outer" electron is ionized first may give a significant 
contribution only when E is of order of the "outer" electron binding energy, i.e. I\. From the 
other hand, it is obvious that the shake-off mechanism may give a significant contribution 
only in the case when the first ionized electron velocity is much larger than the velocity 
distribution width for the second one. Indeed, when the velocity of the first ionized electron 
is comparable to the velocity of the bound one, then the potential affecting to the second 
electron changes slowly, and the second electron adiabatically turns to the bound state with 
the same symmetry as its initial state due to the well-known theorem, therefore its ionization 
can not proceed. Therefore, when the "outer" electron is ionized at first then the shake-off 
mechanism is impossible. But the possibility of the final state scattering is also extremely 
small if the "outer" electron is ionized at first because it can not have an energy enough for 
the ejection of the strongly bounded "inner" electron. Consequently, the double ionization is 
possible if only the "inner" electron is ionized at first. In this case the shake-off mechanism 
is significant only at emission energies much larger than the first ionization potential. From 
the comparing of 7 for the double photoionization with 7 for (e, 2e) in [8] it is clear that 
the angular interelectron correlations are stronger for the shake-off process then for the final 
state scattering. Summarizing all written above we can conclude that 7 decreases with 
the energy decreasing for both processes taken individually, and the 7 increasing at low 
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FIG. 8: The 2DCS as a function of the interelectron angle #12 for the photoionization of He(3s 
for a) E=l eV; b) E=5 eV; and c) E=20 eV: "exact" TDS results (solid line), Gaussian fitting 
cr( 2 ) (dashed line) and the Gaussian fitting of |a 5 | 2 (dotted line). 
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energies in the Figs. [3j [6] and |8] results from the "switching off" of the shake-off process. 
The distinction from the He in the ground state is just the sharpness of this switching off 
because of the narrow momentum spectrum of the "outer" electron. It may be predicted, 
that in the double photoionization of targets with Ji much less than I 2 (like alkali metal 
atoms), j(E) should behave similarly the one for H~ on Fig. [3J 

V. CONCLUSION 

We calculated ab initio the Gaussian width 7 dependence on the electrons energy E in 
the double photoionization of the negative Hydrogen ion H~, the Helium in the ground Is 2 
state and 2s *S and 3s 1 S excited states. For the He(ls 2 ) photoionization our results are in 
the excellent agreement with experimental data but indicate that the well-known Wannier 
threshold law 7 oc E 1 ^ is not correct at energies attainable in modern experiments. It is 
shown that the 7 dependence on the energy is much better described by the curve obtained 
by Kasansky and Ostrovsky when the electrons interaction with the nucleus was taken into 
cosideration. We have shown that for all considered targets with the strongly asymmetrical 
initial state configuration j(E) has the region of the decreasing when E has the order of the 
magnitude of the first ionization potential I\. We suppose that this effect is connected with 
the rapid alternation from the dominating of the shake-off mechanism to the less-correlated 
final state scattering. Also, it was demonstrated that the cross section dependence on the 
interelectron energy for these targets is strongly non-Gaussian at low energies and has a 
number of the peaks equal to the number of the initial state radial nodes, and it reveals the 
new abilities for the targets electron structure qualitative analysis. 
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